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A holographic description on antiferromagnetic quantum phase transition (QPT) induced by mag- 
netic field and the criticality in the vicinity of quantum critical point (QCP) have been investigated 
numerically recently. In this paper, we show that the properties of QPT in this holographic model 
are governed by a CFT dual to the emergent AdSz2 in the IR region, which confirms that the dual 
boundary theory is a strong coupling theory with dynamic exponent z = 2 and logarithmic correc- 
tions appear. We also compare them with the results from Hertz model by solving RG equation 
at its upper critical dimension and with some experimental data from pyrochlores Er2—22 Y22TizO7 


and BiCoPOs. 


I. INTRODUCTION 


Quantum phase transitions (QPTs) happen at zero 
temperature. Such a phase transition and the critical be- 
havior in the vicinity of the corresponding quantum crit- 
ical points (QCPs) have attracted a great deal of interest 
both in theory and experiment sides recently [1—4]. In 
phase diagram, A QCP separates an ordered phase from 
a disordered phase at zero temperature. In contrast to 
classical phase transition at T > 0 where thermal fluctua- 
tions play important role, QPTs are induced by quantum 
fluctuations associated with the Heisenberg uncertainty 
and driven by a control parameter in Hamiltonian rather 
than temperature. Usually, the control parameter could 
be composition, magnetic field, or pressure, etc.. In con- 
densed matter physics, such a quantum criticality is con- 
sidered to be an important mechanism for some of the 
most interesting phenomena, especially in itinerant elec- 
tronic systems [5, 6] and other phenomenons involving 
strongly correlated many-body systems [7, 8]. However, 
the complete theoretical descriptions valid at all the en- 
ergy (or temperature) region are still lacking. 

One of the most interesting and intensively dis- 
cussed QPTs is the ordered-disordered QPT in antifer- 
romagnetic materials induced by magnetic field, such 
as Cu2(C5H12N2)2Cl4 [9], KCuCl; [10], TlCuCls [11], 
BiCoPOs [12] and Er2—-2xY2xTi207 [13]. A schematic 
phase diagram in the vicinity of a QPT as a function 
of magnetic field is shown in Fig 1, where we plot the 
Néel temperature Ty(B) when B < B, and the spin 
gap A(B) at zero temperature when B > Be. In the 
antiferromagnetic (AF) phase (ordered phase), the uni- 
form spontaneously staggered magnetization is accom- 
panied by a long-ranged spin ordering, which persists 
up to a finite transition temperature Ty. When mag- 
netic field increases, the transition temperature Ty is 
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suppressed. When the magnetic field arrives at its crit- 
ical value, B = Be, the transition temperature is sup- 
pressed to zero. When B > Be and T = 0, the quantum 
disordered (QD) phase has gapped magnetic excitations. 
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FIG. 1. Schematic phase diagram for a QPT, showing the 
Néel temperature and spin gap as functions of magnetic 
field [14]. For a magnetic QPT, the characteristic energy 
scales in the QD (quantum disordered) and AF phases are, 
respectively, the spin gap A and Néel temperature Ty, both 
vanish at the QCP. 


Another interesting feature in QPTs involving strong 
coupling is so called “hyperscaling violation” when the ef- 
fective dimension is larger than or equal to de, the upper 
critical dimension of corresponding field theory. In gen- 
eral for the case of dimension less than de, in the vicinity 
of QCP, there are the hyperscaling relations for free en- 
ergy density F(b, T), Néel temperature Ty, correlation 
length € and characteristic energy A with respect to the 
tuning parameter b = B/B,— 1 as follows, 


F,(6,T) =") F (bA, T*), b> 0 (1a) 
Ty x (—b)”, b30-, (1b) 
Each’, Axb” 6307, (1c) 


where d is the spatial dimension of the system, z is 


called dynamic exponent, ~ and v are two positive crit- 
ical exponents depending on model. F,(b, T) is the crit- 
ical contribution to the free energy density defined by 
F.(b,T) = F(b, T) — Freg(b, T) and A is an arbitrary scale 
factor. The effective dimension of the system at QPT 
then is def = d+ z. The scaling relations (1c) are al- 
ways valid. But the naive scaling relation (la) is valid 
only when the effective dimension is less than the upper 
critical dimension, i.e., deg < de [6]. In particular, when 
d = z = 2, the hyperscaling relations (1a) and (1b) are 
no longer valid. In that case, the quantum critical theory 
is not in a weak-coupling region in general. So for this 
case, a description of a strongly coupled effective field 
theory is called for near the QPTs [6]. 


In order to study and characterize strongly coupled 
quantum critical systems, some new methods are needed. 
The gauge/gravity duality or AdS/CFT correspondence 
provides us with a promising approach [15-18]. This du- 
ality relates a weak coupling gravitational theory in a 
(d + 1)-dimensional asymptotically anti de-Sitter (AdS) 
space-time to a d-dimensional strong coupling conformal 
field theory (CFT) in the AdS boundary. In recent years, 
this duality has been extensively applied into condensed 
matter systems. Some remarkable progresses have been 
made in this direction, including holographic superfluids 
(superconductors) [19, 20] (for a recent review, see [21]), 
(non-)Fermi liquids [22-24] and so on. The models in 
the AdS/CFT frame for ferromagnetism/paramagnetism 
and anti-ferromagnetism/paramagnetism phase transi- 
tions have also been proposed in [25, 26]. Especially 
in Ref. [26], we found that the antiferroamgnetic tran- 
sition temperature Ty is indeed suppressed by an exter- 
nal magnetic field and tends to zero when the magnetic 
field approaches to a critical value Be. In this way we 
realized a holographic description of antiferromagnetic 
quantum phase transition induced by an external mag- 
netic field. Since the model involves rank-2 fields, the 
problems about ghost and causal violation may appear, 
which have not been discussed. Very recently, a modified 
model based on the previous works has been proposed in 
Ref. [27], which is shown to be ghost free and causal viola- 
tion is absent, but keeps all the significant results in the 
previous works qualitatively. For this modified model, 
a numerical investigation about antiferromagnetic quan- 
tum phase transition induced by an external magnetic 
field was first presented in Ref. [28]. 


In this paper, we will further elaborate on the magnetic 
induced QPT and study its critical behavior in some de- 
tails. In section II, we will first review some basic prop- 
erties of holographic antiferromagnetic model by combin- 
ing the ideas in Refs. [26, 27] and study the behavior of 
Néel temperature Ty with respect to the external mag- 
netic field, which shows that there is a critical magnetic 
field where quantum criticality appears and can be fit- 
ted well by adding a logarithmic correction to the usual 
power law form. In section III, we will show there is an 
emergent AdSə scaling limit which governs the low fre- 
quency behaviors in the vicinity of QCP. By this emer- 


gent AdS2 geometry, we confirm the numerical results in 
section II and Ref. [28]. In section IV, we give a dual field 
computation from perturbation RG equation inspired by 
holographic results, which gives a confirmation about our 
holographic results. We will also discuss what our model 
can tell us about some real materials in this section. 


II. ANTIFERROMAGNETIC MODEL AND 
CRITICAL TEMPERATURE 


Antiferromagnetic phase in materials is a magnetic or- 
dered phase without net magnetization. For the simplest 
case, there are two different sub-lattices which have two 
spontaneous magnetic moments with the same magni- 
tude but opposite direction. The order parameter then is 
the staggered magnetization, i.e., the difference between 
two different magnetic moments in two sub-lattices. The 
same as in the ferromagnetic phase, the time reversal 
transformation is also broken spontaneously in the anti- 
ferromagnetic phase. 


A. Holographic antiferromagnetic model 


In Ref. [25], we introduced an antisymmetric tensor 
field coupling with U(1) strength field to describe the 
phase which has one uniform spontaneous magnetization. 
In order to describe a spontaneous magnetic ordered 
phase which has a staggered magnetization, we introduce 
two antisymmetric tenser fields. Following Ref. [26] and 
Ref. [27], the action we will consider here takes the form 

1 


S = 


6 V 
~ OK? dtr —glR+5—F* Fu +A°(L1+L2+L12)], 


(2) 


where 


k 
Lz = EMOW MO, 


1 
La) =F 


1 
+ sO Fy + JV(M™), 


2 
(a)\ nvr (qa) M yy (anv grla) 
(aM) (dM) yg + 7M w MG, (3) 


uv 
VM) = CU p MOH), a=1,2. 


L is the radius of AdS space, 2k? = 167G with G the 
Newton constant. * is the Hodge star dual operator. In 
this model, k, m? and J are all model parameters with 
J < 0. à? characterizes the back reaction of the two 
polarization fields M, (a) with a = 1,2 to the background 
geometry, and Lj describes the interaction between two 
polarization fields. The equations of motion for polariza- 
tion fields read 
3V7V MO) — mM — kM® — IV'(a), — Fun = 0. 
(4) 
Here (a,b) = (1,2) or (2,1) and Va), = 
(*M@,,M(77)*M(),,,. In the probe limit of A > 0, 


V 


we can neglect the back reaction of the two polarization 
fields on the background geometry. The background we 
will consider is a dyonic Reissner-Nordstr6m-AdS black 
brane solution of the Einstein-Maxwell theory with a neg- 
ative cosmological constant, and the metric reads [29] 


ds? = r?(—f(r)dt? + dx? + dy?) sui 
rfr) (5) 
1 + py? + B2 ww + B? 
fo) =1- * 


Here both the black brane horizon and AdS radius have 
been set to be unity. The temperature of the black brane 
is 


T = —(3-?— B?), (6) 


For the solution (5), the corresponding gauge potential is 
A, = u(1— 1/r)dt+ Bady. Here p is the chemical poten- 
tial and B can be viewed as the external magnetic field 
in the dual boundary field theory. The zero temperature 
limit then corresponds to B? + u? = 3. 

In order to describe antiferromagnetic phase transition, 
following Ref. [26], we consider M{® , M£) (a = 1,2) and 
define 


1 1 2 
a= MG + MG )), B==(MY) —M&), 


pi = (M MY + mM), p= (Me MÊ). 


N 


1 
2 
Then different values of œ and 8 correspond to different 


magnetic phases. The staggered magnetization NÝ and 
total magnetization N can be defined as, 


Ni/wW=- 
(8) 
N/X =- 


When external magnetic field B = 0, the antiferromag- 
netic phase corresponds to the phase with nonzero stag- 
gered magnetization but zero total magnetization den- 
sity. 

Put the expressions (7) into equations (4), we have the 
equations for a and £ as 


” pa AJ (pt + p3) m? +k 
a+ - [Ae py Jo 
8JpipeB B 
ep F PCr =0, 
IRI 2 2 2 (9) 
BU + fe a AJ (pt + p>) 4 m* —k B 
F PI TJ 
4 Su Pipes = 


r? f 


with 
— r?[(m? — k)r* — 45 (a? + B?)|u 
Pi = 16.F2(a2 — B2)? — 8J (a? + B2)m2rt + (mt — kyr" 
E 8Juapbr? 
Po = 16J2(a2 — B2)2 — 8J (a? + B2)\m2rt + (m4 — kyr 
(10) 


The behavior of the solutions of Eqs. (9) in the UV 
region (near the AdS boundary) depends on the value of 
m? +k. When m? +k = 0, the asymptotic solutions will 
have a logarithmic term, we will not consider this case in 
present paper. On the other hand, when m? + k 4 0, we 
have the solution: 


B 
— „UV ,.(1+61)/2 UV pir) aM 
auy =Q} r FAL T mR 
Buv = pUV rl (1+62)/2 + BUY (1-62)/2 (11) 
ôi = V1 +4k+4m?, 62 = V1—4k+4m?, 
where a¥" and BY" are all finite constants. We require 


that the condensation happens spontaneously when B = 
0, by the inspirit of AdS/CFT, which gives boundary 
conditions such that 


a =p = 6. (12) 


B. Conditions for antiferromagnetic phase 
transition when B=0 


In our model, there are three parameters J,m? and 


k. In order that the dual boundary theory can describe 
antiferromagnetic system, these parameters need satisfy 
some conditions. 

First, from the asymptotic behaviors in Eqs. (11), w 
see that the BF-bound at the boundary requires that m? 
and k satisfy 1+4m? > 4|k|. On the other hand, because 
of J < 0, in order to make pı and pə finite in the region 
(rp, 00), we need following conditions, 


m? > |k| å > 1, and ô > 1. (13) 


If conditions (13) are broken, we can see that the compo- 
nents pı, p2 will diverge at somewhere, which leads that 
there is some point where the energy-momentum tensor 
associated with these tensor fields diverges. So in this 
case, the bulk geometry becomes unstable and the probe 
limit is broken. 

We can explain the conditions (12) and conditions (13) 
in other way. By the definitions (7), we see that the terms 
of a2’, BYY and B all appear in the asymptotic solu- 


tions of MS) and MS). Since we need that the system 
is dominated by external magnetic field when B Æ 0, 
so the magnetic term should be the leading term near 
the boundary. This can be reached only when both the 
boundary conditions (12) and conditions (13) are satis- 
fied. 


Second, as pointed out in Ref. [26], in order to get an 
antiferromagnetic phase transition below a critical tem- 
perature when B = 0, we should achieve that 6 can con- 
dense spontaneously, but œ cannot, or in other words, 
the system is stable for a but not for 6. In RN-AdS 
black brane background, this can be achieved if the AdS2 
BF-bound is satisfied for œa but violated for 6 near the 
horizon when T = 0. As in the critical cases for instabil- 
ity, we can treat a and 8 both infinitesimal, so we get a 
linearized equations for them, 


gente... — taso, 


f 2f TJ 
(14) 
IRI 2 2 


with pı = u/|(m? + k)r?]. Then zero temperature corre- 
sponds to p = V3. Let rą = 1/(r—1), when rą — 00, we 
have following asymptotic solutions for a and £, 


I —(14+67)/2 f —(1-67)/2 
arr = of Bry! v/ + of By Z' 1) ; 

TR,,—(1+63)/2 IR —(1—ô3)/2 
B = Tx ( 2)/ B Tx 2)/ ; 


: ; (15) 
ï =4/1+2(k +m? + 4Jpiy)/3, 
> =4/1 +2(m? — k + 4Jp?o)/3. 


Here pio = V3/(m? + k), alF and 61? are all finite 
constants. As a result, to have the antiferromagnetic 
instability requires that the parameters satisfy following 
conditions, 


[ 1 +2(k +m? + 4Jp?,)/3 > 0, (16) 
1+ 2(m? — k +4Jp?,)/3 <0.’ 
which gives following constraint, 

Jè (k, m?) < J < JZ hen), and k > U. (17) 


with JE (k, m?) = —(m? + k)? (m? + 3/2 + k)/12. 

When external magnetic field B = 0, for given mass 
square m7, the zero temperature phase of the dual bound- 
ary then depends on the parameters J and k. A typical 
example for m? = 1 is plotted in Fig. 2. The red region 
corresponds to the case that parameters satisfy the con- 
ditions (13) and (17), which is what we will consider in 
this paper. 


C. Néel temperature and quantum criticality 


When the magnetic field is absent, i.e., B = 0, under 
the conidtions (12), (13) and (17), the solution of a is 
always zero, which means there does not exist sponta- 
neous magnetization. But there is a critical temperature 
Tyo, lower than which a nonvanishing Ø begins to ap- 
pear, which shows the system transits into antiferromag- 
netic phase and gives a nonzero spontaneous staggered 
magnetization. 
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FIG. 2. Zero temperature phase diagram when B = 0 and 
m? = 1. Here FM, AFM, PM and Ferri stand for ferromag- 
netic phase, antiferromagnetic phase, paramagnetic phase and 
ferrimagnetic phase, respectively. 


On the other hand, if turn on an external magnetic 
field B, we have a # 0. When the external magnetic 
field is small, the value of |a| will increase with the in- 
creasing of B. Since the equation for 8 couples with a, 
the staggered magnetization is also influenced by mag- 
netic field. This influence can be seen clearly from the 
effective mass square of 8 in the IR region, namely, the 
near horizon region. Near the horizon, we have, 


4Jp? 2 


4J? u22 
ee ee, 64J ura 


(e+ m2 (m2 — by 

(18) 
Under the restriction (13), we can see that, for small mag- 
netic field, the effective mass square of 3 in the IR region 
will be increased by magnetic field B, which leads the 
nonzero solution of 8 to appear in lower and lower tem- 
perature with increasing the magnetic field. Thus the 
antiferromagnetic critical temperature Ty will be sup- 
pressed by external magnetic field. 


=hF 


2 
M Beff 


When the external magnetic is not very small, because 
of the nonlinear coupling between 8 and a, it is not easy 
to directly give the relation between effective mass square 
of 8 and the external magnetic field. Instead, the behav- 
ior of Ty with respect to magnetic field can be found nu- 
merically. Near the critical temperature, the staggered 
magnetization is very small, i.e., 6 is a small quantity. In 
that case we can neglect the nonlinear terms of 6 in the 
equations (9), which leads to 


" ai Ja Moff B 

I r o 

ral 2 (19) 
B" + FB Mpe g = 


with 


mog = 4Jpi +m? +k, 

Meek =m? — k + 4Jp? + 8Jpipro’, 

r? y3- Bia w) 
(m? + k)r4 — 4Ja?’ 

m 8Jr? v3 — B2 — 4nT 


R= (4Ja? — m?r4)? — k?r8 ` 


Pı = 


When T Æ 0, in the IR region, the regularity of the 
solution gives following initial conditions 


a = OM? er E B 
ArT ’ 


2 
BM beer 


td 
P= ar 


(21) 


at the horizon. Without loss of the generality, we can set 
B(rn) = 1 due to the linearity of the equation of £. 
When Ty = 0, the initial conditions (21) imply that 


B= am oglr, and Meeetlrn =0, (22) 


at the horizon. For given parameters J,m? and k, these 
two equations give the solution a(r;,) = a. and B = Be. 
Here az is the initial value of a at the IR fixed point and 
Be is the critical magnetic field where the QPT occurs. 

In Fig. 3, we plot the Néel temperature Ty in the 
case of B < Be. Because different parameters satis- 
fying restrictions (13) and (17) give qualitatively same 
results, we here just show a typical example by taking 
m? = 1,k = 7/8 and J = —2/3. With the increasing 
the magnetic field from zero to Be, the Néel temperature 
decreases from Tyo to zero. For small B, numerical re- 
sults show that Ty — Tyo x B?. When magnetic field 
approaches to Be, we find that the Néel temperature is 
fitted well by following relation 


Ty /InTy ~ —0.4074(1 — B/B.). (23) 


Here Ty = Ty /Tno and Tyo is the Néel temperature in 
the case of B = 0. 

The numerical results shown in Fig.3 and the rela- 
tion (23) show that there is indeed a critical magnetic 
field Be, at which the critical temperature for staggered 
magnetization, i.e. 3, is zero. So a quantum phase tran- 
sition occurs, which separates antiferromagnetic phase 
and a quantum disordered phase at zero temperature. 
We can also obtain other critical properties in the vicin- 
ity of QCP by numerical methods [28]. Especially, in the 
vicinity of QCP, this model gives a dispersion relation for 
quasi-particle such as w ~ q? with energy w and moment 
q, which shows that the boundary critical theory is a 
strong coupling theory with dynamic exponent z = 2 [6]. 
This agrees with the results from condensed matter the- 
ory about antiferromagnetic metal [6] and from the holo- 
graphic model proposed by Ref. [30]. 
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FIG. 3. The antiferromagnetic critical temperature Ty versus 
the magnetic field B. (a) In the whole region of 0 < B < Be. 
(b) In the region of B < Be. (c) In the region of 1— B/B. > 
oO". 


The behavior of the Néel temperature in the vicinity 
of QCP is quite non-trivial and needs to be further un- 
derstood. The relation (23) is not the usual power-law 
behavior in (1b) or square-root form. However, it agrees 
with the predication in the 2-D QPTs in strong coupling 
case [2, 6, 31]. The d = z = 2 quantum critical theory is 
in general not in a weak coupling region for any T > 0 
and a strongly coupled effective classical model emerges 
that can be used to determine the critical dynamics [32]. 
As pointed out in Ref. [28], it is just the reason that 
the AdS/CFT correspondence is very suitable for this 
description. In addition, the original numerical results 
in Ref. [28] need to be confirmed, especial for the log- 
arithmic correction scaling rule Eq. (23). And we also 
need some more careful investigation to provide insight 
into physical properties in the vicinity of QCP in this 
holographic model. The most important is to clarify the 
relationship between this holographic model and the tra- 
ditional theories of QPT and to know what we can learn 
about real materials in experiment from our holographic 
results. These are our goals in the following sections. 


Ill. EMERGENT IR GEOMETRY AND THE 
BEHAVIORS NEAR QCP 


In order to understand the scaling relation (23) and 
other numerical results presented in Ref. [28], we now 
pay attention to the background geometry in the zero 
temperature case. As the similar situation in Ref. [33], 
we will see that our numerical results about the behav- 
iors in the vicinity of QCP are just controlled by the IR 
geometry which is governed by an emergent AdS geom- 
etry. 


A. IR fixed point and critical magnetic field 


In the vicinity of QCP, the system is dominated by the 
features of IR region, where an AdS geometry emerges. 
Following Ref. [33], we introduce a new length scale 7, 


R+B? = 34 (24) 


Then the temperature of the black brane background can 
be written as, 
pa (1—7*) (25) 
= —(1 -f 
4T 
Using this new scale, we define a new radial coordinate 
€ and time 7 as, 
~ A ~ or 
r-Tr=—, 1-r=— 
6E 6&0 
Use the variable €, we can define an inner IR boundary at 
E — ĉo and an outer IR boundary at € — 0* (see Fig. 4). 
The line element can be written as, 


, t=d7'7, A340. (26) 


1 2 E? 
ds? = — |—(1— =,)dr? + (1 — 35) 1 de" | + (dr? +dy? 
(27) 
Then the radial function f(r) is, 
2 E2 
= — (1 - >). 28 
f= l-h) (28) 
And the temperature with respect to T is, 
~ 1 
Pa) = l 
À A (29) 


Here the case of zero temperature corresponds to an in- 
finity €). Note that in the scaling limit (26), finite + 
corresponds to the long time limit of the original time 
coordinate. Thus in the language of the boundary field 
theory, the solution (27) corresponds to the low frequency 
limit. 

Since the radial coordinate can be considered as the 
resolution scale for the dual theory. The evolution of the 
geometry and the fields propagating therein along this 
radial direction represent the RG-flow of the dual field 
theory. Using the coordinate transformations (26), we 
can rewrite the equations (9) in the limit of A — 0 as 


da? = X% da E(B+am2g) , 


a Q-PE oPR =E) 

dp? E 2g dp = Eo MBerb _9 

dẹ — E72 dE 6ER(EG-—E?) 
In the zero temperature case, which corresponds to éo > 
oo, Eqs. (30) become 


(30) 


da? B+ Mca es 


(31) 
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FIG. 4. A schematic figure for the inner IR region and outer 
bulk region. The horizon is located at € = £o. The near hori- 
zon region is o < € < 0, where an AdS2 geometry emerges. 
The bulk region is not covered by €-coordinate. The solutions 
in the IR region and the bulk region should be matched at 
€=0. 


With the boundary conditions that œ and ĝ are both 
finite at inner IR boundary and outer IR boundary or 
the existence of IR fixed points for a and p in the IR 
region, these two equations have regular solutions only 
when a and 8 are both constants. As a result, we have, 


B+ moga =0, Mog = 0. (32) 


This is just what we have obtained in Eqs. (22). Note 
that in the case, a, 3,p, and p2 are all constants in IR 
region, SO aM*,.¢ and Mog are also constants in IR re- 
gion. With the restrictions that m? > |k| and J satisfies 
Eqs. (17), we can conclude that when B < Be, the solu- 
tion with 8 Æ 0 exists. But when B > Be, there is only a 
trivial solution with 8 = 0, by matching into bulk region, 
which will give a trivial solution such that (r) = 0 in 
the whole bulk region. Thus there is a phase transition 
at B = B, at zero temperature, which divides an AFM 
phase from a quantum disordered phase. 


B. Linear perturbations and dynamic exponent 


In order to compute the magnetic fluctuations in the 
vicinity of QCP. We need to turn on the perturbations of 
two polarization fields. Because of the non-linear term in 
equations (4), all the components couple with each other. 
As a result, we need consider the perturbations for all the 
components, 1.€., 


SMI) = Oeit), (u,v) A (r, y), (t,x) 
5M) = iC eH), (u,v) = (r, y), (t, 2). 


Substituting the ansatz (33) into (4) and compute up to 
the linear order of €, we can get the equations for per- 
turbations. When w,q 4 0, however, we can’t decouple 
all the equations. In general, therefore we needs to solve 
the 6 x 2 components together in order to determine the 
dispersion relation. But if we consider the behavior for 
small frequency and wave vector in the vicinity of QCP, 


i.e., w + 0 and q > 0, the problem can be simplified. In 


that case, the equations for c® and Cry (a) decouple from 
the others and can be neglected. Then s we obtain 4 x 2 
coupled equations, 

[(a? + mr? - aM? CY + FC 


P + kro = 0 


2 470 MS ME 
2 (a) (b) — 
aeo oO SE ee 
(nt Gap) 8 + T 
a b a a a 
ay” mC kO . AIC MSM _ A 
ty r2 f rf r4 f a 
(a)! 2 (a)2 
oo” i S Cri a wr m? + 4JM,.. co’ 
f r2f2 a y 
a b (a a 
Cga pOL BIOP MT ME _ 4 
rif? r2f Me i 
(34) 
with (a,b) = (1,2) or (2,1). In the paramagnetic phase, 


ie, MẸ) = MÊ and M® = MÊ, the second and 
the third equations in (34) show that Of = = O a nd 


cH = 0. Let 8 = (cQ — Cc) /2, then we obtain, 
xy P'E |w-QPf Moen 
1 
B P + rr ee B+O(q 4) = 0. (35) 
with 
64J2 A na 
Ga yi (36) 


(m? — k?)r4 — 4Ja?]? 

As q is infinitesimal, we can neglect the high order term 
O(q*). Note that no matter how small w and q are, we 
here can’t neglect the (w? — Qq?)/(r?f?) term, because 
it diverges when f = 0. 

In order to analyze the antiferromagnetic dispersion 
relation near the QCP for small frequency w and wave 
vector q, we can take the IR-UV matching method pro- 
posed in Ref. [33]. We first compute the retarded Green’s 
function in the IR region. Use the metric (27) in the limit 
of ĉo — co and scale frequency as w = Aw, we obtain, 


PP 


Geet @ - od? /662)8 = (37) 


where we have used the results (32). Qo is the value of 
Q in the IR region. Because a is a constant in IR region, 
Q is also a constant in IR region. Solving this equation 
with the ingoing condition at € — ov, we have, 


Bx JeH Ge). (38) 


Here H, P(x) is the first Hankel function with order v = 
\/9 + 6Qoq?/6. Using the properties of Hankel function 


and the fact that q < 1, we can find that the asymptotic 
solution for 6 in the outer IR region is, 


B x 67 00/2 4 G(@)g1+ 209" /2 (39) 
with the IR retarded Green’s function G(@) œx 0209/3, 


Then using the method in Ref. [33], one can find that the 
retarded Green’s function can be expanded in small w as, 


G(w) x 
DO + wo) + O(w?) + G(w)(bO + wo + O(w?)) 
x + wal) + O(w?) + G(w) (a twa + O(w?)) 


(40) 


with b and a\? are the analytical functions of q?. By ne- 
glecting the irrelevant parts in (40), the retarded Green’s 
function can be written as the following form for small w 
and q, 


Z 


CC E. ER 
( ) —¢q? +ayw+ anu oe /3 


(41) 


with some constant Z. Considering the fact that when 
q = 0, the retarded Green’s function has a pole at w = 0, 
we have ag = 0. Then the dispersion relation is deter- 
mined by the equation, 

— Qog? + aw = 0, (42) 
which tells us w ~ q?. Thus we obtain the dynamic 
exponent z = 2. 

One should note that, in general, the coefficients ap- 
pearing in (40) have scaling dimensions, and so dose for 
a, in (42). From the expression (42), we find that a, has 
scaling dimension 1. Thus under the scaling transfor- 
mation such that r > Ar, (t,2,y) 3 A7l(t,2,y), which 
induces the transformations such that w > Aw,q > Aq 
and a; — Aaj, the equation (42) is invariant. As a result, 
the dynamic exponent z = 2 is compatible with scaling 
symmetry. 


C. The scaling relation of Ty with B near the QCP 


Now let us find the the scaling relation of Ty with B 
in the region of B + Bz. The method we will take is 
similar to the one as in the previous subsection. We first 
compute the Green’s function at finite temperature in 
the IR region, and then match it with the bulk solution. 

In the IR region, near the fixed point B = Be + ôB 
with 6B — 07, we can neglect the non-linear term of 3 
and can get the equations for finite £o, 


da? 2 da 

d BO de 
dp? 2g dp 
dẹ? @-Pd 


[B+ Mera] 
6E E) 

_ ESM oon 

6A =e) 


(43) 


Under the boundary conditions that œ and ĝ are both 
finite at € = 0 and ĉo, the solutions for Eqs. (43) are 
constants. Then we have solution such that, 

B.+6B+m .~«a=0, 6=0, (44) 
Using Eqs. (22), up to the order of ôB, we find that 
Eqs. (44) give solution a = ao = a, + ĝa with 


ôB 
ĉa = — — ; 45 
(mef + Om? o¢/OQ) l ) 


aA=Ac 
In order to compute the retarded Green’s function, we 
consider the deviations from the fixed point with 6 = 
0+ 8, and a = ao +a ,. Note that 6B — 07, up to the 
linear order of 6, and a,, we have the equations of the 
fluctuations 

da? 2E day 


& (men + Om? .¢/Oa)ay 
bet — £2) i 
doè dbr _ mher _ sapı __ 
dé? ge — E dE ða 6E (ER — E) 
Here effective mass terms and their derivative are taken 


the value at @ = ao. The solutions for a; and 8; can be 
expressed as hypergeometric functions as, 


dE -EdE 
dB? 2 dp 


(46) 


1t5a (a 346, Ltdg Oa 
a =) ? carl 7 og = Sieh 
£146 Jip i43 5 
a= Dey cP ne Ae BS 
l (47) 


with ĝa = y1++2(m2.ş +3mMź g/0a)/3 and dg = 


1+ ZOM? es ða)ða. These hypergeometric functions 


aeff 


have logarithmic divergency when €/& — 1. So the co- 
efficients oe) and of ) should counteract these deter- 


gency. Then we can get the retarded Green’s functions 
in the outer IR region near € —> 0 as, 


3 4 da\p(l 4 bo — Sa 
Gaa(F) = (any a+ Sta + ECB) 
4 


T( 


Now we can use the IR-UV matching method to express 
the UV Green’s functions. Similar to the expression (40), 
we can get the Green’s functions for a; and 81 with some 
coefficients a”) and b{” which should be analytical func- 
tions of 6B. Note that two retarded Green’s functions 
have poles at 6B = T = 0, we have, 


ÖB + ST + Gaa lT)” 
ôB + CPT + Geg (T)c) , 


aa p 


(49) 
Gap 


Here Gaa is irrelevant and can be neglected. One should 
note cP ) are also the functions of of) because the equa- 
tion for 6 in the bulk depends on a but the equation for 
a in the bulk is independent on £8. In the limit of 6B = 0, 
we will find 6g = 1, so the of) term and P term de- 
generate. In that case, as the similar as the case pointed 
out in Ref. [33], a logarithmic term appears. And the 
relevant terms in Green’s functions read 


Zı Z2 


Gaa = ————, Gas = — n 
SBOT 5B 4 PO TIT 


(50) 


with two constants Z; and Z2. The poles for them are 
defined by, 


B+ ®T=0, B+) TInT=0. (51) 


Now note the fact that 6B — ôB, c > -c will 
lead a + —a, which will not affect the solution of 2, 
cf ) must be a function of fo? Then in the case of 
T — 0 and ôB — 0, there is a self-consistent ansatz such 
as of) a dye?” + O(c) for small of? and finite 
constant d;.' Then we can see that the relevant parts of 
the equations (51) can be expressed as, 


5B+0°T =0, 5B +c" dTMT=0. (52) 


Solving these equations, we find ôB x T/InT. Asa 
result we have b = 6B/B. x T/I\nT and analytically 
show the relation (23). 


D. Energy gap and correlation length in the 
quantum disordered phase 


In Ref. [28], numerical results show that there is a 
gapped antiferromagnetic excitation in the quantum dis- 
ordered phase. The energy gap and correlation length 
still obey the power law form (1c) and (1b) with respect 
to the tuning parameter b in the vicinity of QPT. These 
results can also be understood from the point of view of 
the emergent AdS2 geometry. 

In order to compute the energy gap of antiferromag- 
netic excitation, we should turn on a perturbation for 8 
with frequency w. In the region of quantum disordered 
phase with T = 0 and 6B = B — Be > 07, using the 
results in the previous section, the background and per- 
turbation equations for polarization fields read, 


1 In general, the expansion of P at small ee has a constant 


term such that do. This constant term must be zero, otherwise, 
the equations (51) will not have solution if 6B 4 0. 


and 8 = 0. The real frequency w giving the pole for 
Green’s function Ggg(B) leads to the gapped quasi- 
particle excitation. To compute this Green’s function, we 
take the similar method as in the previous subsections. 
In the IR region, considering the AdS2 scaling limit at 
zero temperature, the equation (53) can be written as, 


Pa B,+6B+m2 ga _ ð 


Je oe (54a) 


(54b) 
Up to the linear order of 5B, the equation (54a) gives 

the following regular solution with the IR fixed point, 
a = Qe + ĝa. (55) 


Here ôa is defined as (45). The equation for B can be 
written as 


28 ~2 OM ort ôa in 
we fe -( a æo 0 


c 


Comparing it with the equation (37) in the subsection 
IIIB and note the relationship in (45), one can easily 


find that the Green’s function for 8 has the following 
form, 


Z3 
G = ————_ 57 
La ew (57) 
for small frequency w and 6B — 0+ with two nonzero 
constants e; and Z3. Thus we obtain the relation between 
energy gap A and magnetic field B as 


AxB- Be, (58) 


As to correlation length, it can be computed in a simi- 
lar way by setting w = 0 but turning on the perturbation 
for @ with nonzero moment q. In that case, we still have 
equation (54a) but, the equation (54b) is changed to be 


2B 7 Qog? + Moet _ 


eB a 0. (59) 


For small q and ôB, it gives the Green’s function as, 


Z4 Z4 


G PEE PR E ee 60 
ap w) q? + €26B P FIJE ( ) 
where Z4 is another constant and € is called correlation 


length. We then easily obtain 
Eapen (61) 


which is as the same as equation (1b) with v = 1/2. 

We see that results (61) and (59) obey the universal 
scaling relations (1b) for quantum criticality with z = 2, 
ie, A x |B — B.|*”. This can be viewed as the self- 
consistency check for our model and computations. 


IV. CONCLUSIONS AND DISCUSSIONS 
A. Summarize the holographic model 


We have investigated the quantum phase transition 
(QPT) from antiferromagnetic phase (AF) to quantum 
disordered phase (QD) and the criticality in the vicinity 
of quantum critical pointc(QCP) in the antiferroamgen- 
tic materials induced by magnetic field in a holographic 
model proposed in Ref. [26]. In the model, the Néel 
temperature Ty is suppressed by magnetic field B. We 
proved that there is a critical magnetic field Be, at which 
Ty = 0 and a QPT occurs. In order to analytically 
study the critical behavior in the vicinity of QCP, we 
employed the IR-UV matching methods in Ref. [33], and 
the results show that the analytical method confirms the 
numerical results obtained in Ref. [28]. In particular, we 
found that the dynamic exponent z = 2, which means 
that the boundary critical theory is indeed a strong cou- 
pling theory with effective dimension deg = d+z = 4. As 
a result, the hyperscaling law is violated and logarithmic 
corrections appear near the QCP. An interesting observa- 
tion is that the critical behavior of the QCP is governed 
by the AdS2 geometry emerging in the IR region of the 
background geometry. 


B. Compare the results with RG-flow 


In this subsection, we will review the results from field 
theory and compare them with the results from the holo- 
graphic model. Our discussion follows Ref. [6], the details 
can be found there. A similar discussion as what we will 
do but mainly for the case of d = z = 2 can be found in 
Ref. [35]. 

In the original work by Hertz to quantum criticality [2], 
he considered the coupling between fermions and the low- 
energy bosonic field which condenses at the QCP. By in- 
tegrating the fast fermions, the effective field theory is 
the Landau-Ginsburg-Wilson (LGW) ¢* theory with the 
upper critical dimension de = 4— z, where z is the dynam- 
ical exponent. Near the QCP, we can get the perturbative 
RG [34]. Let u, Ty, A and b are the coupling constant, 
condensed temperature for bosonic field, scaling param- 
eter and tuning parameter (in this paper, it’s defined as 
b = (B — B.)/B), respectively. Then the infinitesimal 
RG transformations are given by [34] 


dT 
mm ve 
An eee nes (62b) 
ina Oe 
d L (4-d-zju-u?fa(T,b). (620) 
mma “i zju — u’ fo(T, b). ê 
with some initial value 7 = Ty,b = bo,u = uo. Here 


fi(7,b) and fo(T, 0) are two functions of T, b defined by 
the details of model. At low T close to the critical point 


0.57 


/ 
RA 


-0.5 + 


FIG. 5. The schematic RG-flow diagram of u-b for d+ z > 4. 
There is a Gaussian fixed point at u = b = 0. 


b = 0, fı(T,b) can be expressed as the power series of 
Tx and f2(7,b) is constant. Thus we can set fı(T,b) = 
Ao + AiTX, and f2(T,6) = f2. The RG equations (62) 
have a Gaussian fixed point at T = u = b = 0, which is 
unstable with respect to tuning parameter b, so a phase 
transition may occur at this point (see Fig. 5). Since the 
RG equations are obtained by perturbations, we assume 
that the initial values satisfy Ty, |bo|,|wol < 1. One of 
remarkable features is that this RG equation has a critical 
dimension at de = d+ z = 4. So, in general, when 
d > 4-— z, we can find that the boundary of ordered 
phase Ty (such as the AF transition temperature, see 
Fig. 1), energy scale A in the quantum disordered phase 
and the correlation length / in the vicinity of QCP obey 
the following scaling relations, 

Ty ~ (=b), Awd, lw br”. (63) 
with Y = z/(d+ z — 2), v = 1/2. However, this simple 
power law will lose its efficacy in the case of d = z = 2, 
where the RG equations (62) need special consideration. 
In fact, logarithmic correction on Ty appears. In order 
to see this, we directly solve 7,u from the RG equa- 
tion (62a), 


uo 


T(A) = Ty, u(A) = are are 


(64) 
Next, integrating the scaling Eq. (62b) for b(A) and using 
Eqs. (62a) and (64), we find 


BA =>? [bo + 20x Ei (2 ma + > /(uosa)] (6) 


Uo J2 


Here F(x) is exponential integral function. For large 
real x, Ey(x) ~ e~*/x. In the quantum critical region, 


using the scaling up to A ~ TUR where 7 ~ 1, we have, 


bo 1 1 


W= pe + Fonts) O TnT" 


(66) 
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Thus we see that the scaling relation depends on the value 
of b(A) as A + œ. However, the concrete form of b(A) 
can not be found from the perturbational RG equations 
directly. Because the feature of strong coupling in the 
case of z = d = 2, the perturbational methods can’t give 
complete information of the system. We see that the 
perturbation method fails in this case to give a complete 
result. If b(A) is a finite constant, then we can find a 
linear relation bọ ~ Ty, which is just the usual power 
law. If b(A\) nA > 0 as A = 1//Tw > œ, then we can 
get bo ~ Ty/InTy. Comparing with the usual power 
law, there is an additional logarithmic correction, just as 
the same as what we have obtained in the holographic 
model. 

We see that, with some additional assumption, the per- 
turbational RG equations give the same results as the 
holographic model. In fact, in the case of d = z = 2, 
the low-energy Fermi liquid leads to long-range order- 
parameter interaction. As a result, the coefficients of the 
high-order interactions in the LGW functional diverge, 
leading to an infinite number marginal operators [31]. 
We see that the perturbational field method in this case 
fails to give out a complete description. However, the 
holographic setup can give self-consistent and complete 
predictions. 


C. Apply to real materials 


Now let’s turn our attention into real materials and see 
what we can learn from this holographic model about the 
real materials. Here we focus on ErəTi207, as a kind of 
Magnetic pyrochlore oxides materials, which has amazing 
properties. At low temperatures there exists an antifer- 
romagnetic order arising with a Néel temperature of 1.2 
K. The compound is complex and there together with the 
nearest neighbor exchange, the Dzyaloshinskii-Moriya, 
dipolar and magnetoelastic interactions do exist and play 
an important role in the formation of the low tempera- 
ture antiferromagnetic state. There is spins ordered on 
the vertex-sharing tetrahedra, in non-coplanar structure. 
There order is formed by quantum disorder, i.e., in a 
very unconventional way [13]. To describe this mate- 
rial the XY-like theoretical model where local spins are 
coupled by near-neighbour antiferromagnetic exchange 
and long-range dipole interactions on the pyrochlore lat- 
tice has been proposed [36]. The key role in the model 
is played by the isotropic exchange between transverse 
components of near-neighbour spins and bond-dependent 
exchange anisotropy. The authors have shown that the 
interplay of these interactions dictates the formation of 
the ground state by the quantum disorder mechanism. 
The inclusion of the exchange anisotropy reproduces well 
the observed second-order magnetic phase transition in 
ErəTi207 at zero field. However, in Ref. [37] it was ar- 
gued that the dipolar interactions makes an important 
contribution into the formation of the second order phase 
transition observed [13]. Probably as argued in [36] the 


exchange anisotropy may overcome the dipolar interac- 
tion at zero magnetic field where the tetrahedral mag- 
netic units have no net magnetic moment in the ground 
state. However with the field when the canted states are 
formed the dipolar-dipolar interaction, which certainly 
exists between large magnetic moments associated with 
free Er?+ ions (J = 15/2, gy = 6/5), may become very 
important and can not be neglected. Therefore, the field 
evolution of the antiferromagnetic state and the role of 
the long-range part of the dipolar interactions in such 
an evolution remained unclear and beyond the scope of 
the discussed models. Here we proposed a holographic 
model which effectively absorbs all existing interactions 
in ErəTi207 in a form of strongly interacting tensor fields 
embedded in AdS space of a higher dimension. The extra 
dimension play a role of the renormalization parameter, 
which separates relevant and irrelevant degrees of free- 
dom in the system. The proposed holographic model 
describes well not only the zero temperature thermody- 
namic phase transition and spin wave excitation spec- 
trum in Er2Ti2O7 but also their evolution in the applied 
magnetic field [28]. The model may uncover an intrigu- 
ing picture of many possible coexisting phases, which 
can arise at such complexity of interactions existing in 
ErgTi2gO7 [38]. 


In particular, in Ref. [38] it was found that with mag- 
netic order there arise well defined spin wave excitations. 
With applied magnetic field their spectrum is softening 
(see the spin wave spectrum measured at the field 1.5 
T on the Fig.4 in Ref. [38]). At the studies of an ar- 
ray of magnetic particles where the dipolar interaction is 
usually relevant it was shown that its long range char- 
acter may lead to a formation of a landscape of low- 
energy states [39, 40]. With increasing field the spins 
form canted state and the energy landscape is flattering 
and there may arise continuous switching between these 
states. Therefore such evolution may lead to a quantum 
critical behavior and an appearance of quantum critical 
point QCP. Indeed the analysis of the low temperature 
phase diagram of ErəTi2O7 made in Ref. [38] supports 
an existence of the rich variety of elementary excitations 
as the ground state is tuned by external magnetic field. 
The authors of Ref. [38] also suggested that most dra- 
matic variations of the spectra with field arising at 1.5 T 
may possibly correspond to the formation of a continuous 
QCP. Indeed in the framework of our holographic model 
similar picture emerges. The holographic model predicts 
that the magnetic order continuously changes with the 
field and the QCP arises at critical field Be ~ 1.45 T. 
Our model provides not only qualitative but also quan- 
titative description of the experimental data [12, 13] of 
the Néel temperature lowering with magnetic field (see, 
Fig. 6) and predict the QCP at 1.45T, which is very 
closed to the experimental results. There magnetic spins 
become partially aligned in the direction of the mag- 
netic field. Therefore the system requires less thermal 
energy to destroy the remaining magnetic spins order. 
We found that similar behavior exist in other pyrochlore 
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FIG. 6. The antiferromagnetic critical temperature Ty versus 
the external magnetic field B calculated at the best fitting 
model parameters m? = 1,k = 7/8 and J = —0.71 (solid 
curve) and compared with experimental data for pyrochlore 
compounds: Er2—2sY2sTi207. Black circles corresponds to 
zero doping, x = 0; magenta triangulars - to x = 0.031; green 
crosses - to x = 0.085. The experimental data are from [13] 
and the figure is from [28]. 


compounds which may be obtained from ErəTi2O7 by 
Y+ doping. Both doping and the increasing of magnetic 
field decrease the Néel temperature of the Y*+-doped 
and undoped Er2Ti207 very similar until the QCP arises 
(see, Fig. 6). We see from this figure that the critical tem- 
perature versus critical field curve is independent of the 
diamagnetic doping at the magnetic site. This fact that 
the influence of dilution in (Er;_,Y,)2Ti2O7 solid solu- 
tions is similar to the increase of magnetic field for pure 
Er2Ti2O7 one more time stresses the importance and uni- 
versality of complex magnetic interactions in these sys- 
tems which can be well described in the framework of the 
proposed holographic model. 


Here we also found that the high-field ground state 
does not correspond to a simple disordered paramagnet. 
The high-field state is characterized by the intense dis- 
persive excitation, which dispersion increases with the 
field. This was exactly observed in Ref. [38] where it was 
found that the critical fluctuations of the low-field state 
are accompanied by a precursor of the intense dispersive 
excitation noticed at high field. In Fig. 7, we show the 
magnetic gap energy of Erg_9,Yo,Ti2O7 as a function 
of magnetic dilution (x) and magnetic field [13]. One 
should note that, in Fig. 4 of Ref. [38], the excitation 
gap energy is zero when B < Be and T — 0. But by the 
magnetic heat capacity method in Ref. [13], there is a 
nonzero gap energy (œ 0.18meV) when B < Be (see the 
case of B = 1T in Fig. 7). One can see that this nonzero 
gap energy is nearly independent on the doping x which 
influences the antiferromagnetic electronic properties. So 
this nonzero gap in fact is from the background magnetic 
heat capacity which has nothing to do with the gap en- 
ergy of the antiferromagnetic exaction. After removing 
this background contribution, we can find that the en- 
ergy gap just satisfies the linear relationship (58), which 


FIG. 7. The gap energy,A, of Er2—2s Yor, Ti2zO7 as a function 
of magnetic dilution (x) and magnetic field. The data is from 
Ref. [13]. 


is equivalent to 
A x B/Be-1 (67) 


with A = A/Ty. In the subfigure of Fig. 7 we plot 
the experimental data and the fitting curve by Eq. (67). 
The slope of the fitting curve is 4.2. We can also obtain 
the this slope from the holographic approach under the 
parameters of the best fitting from the Fig. 6, which gives 
that the slope is 5.0. So we see that our holographic 
model can not only fit well the Ty-B behavior but also 
predict the critical magnetic field which is very close to 
the measurement result and the linear behavior of the 
energy gap with magnetic field with the slope nearly to 
the measurements at the same time. 


We also found that holographic model describe well 
magnetic properties of other class of materials, which 
have strong magnetic anisotropy and large magnetic mo- 
ments. The most illustrative example is BiCoPOs; [12, 
28]. The compound consists of the interacting spins 
chains giving the 1D character to the compound. The 
spins have large magnetic moments, about 5.2 upg per 
Co ion measured at 1.5 K [12]. It is much higher than 
the value 3.9 ug expected for Co?* ions having S=3/2. 
Such large value of jeg is rather typical in Co?+ con- 
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taining compounds due to additional orbital magnetic 
moment created in oxygen octahedra existed around the 
Co?* ions. 

Magnetic structure has been determined with the pow- 
der neutron diffraction experiments [41]. It was sug- 
gested there the ferromagnetic pairs of Co?+ ions are 
antiferromagnetically coupled within double chains. Be- 
cause of these large values of the ferromagnetic pairs of 
Co?* ions there exists a dipole-dipole interaction between 
the resulting magnetic moments, which is, probably, sig- 
nificantly, weaker than the exchange interactions. There 
is also strong the spin-orbital interaction, which always 
exists in heavy atoms such as Bi. The competition of all 
these interactions is probably responsible for the devia- 
tion from collinear AF ground state in BiCoPOs [41]. In 
this sense the situation is very similar to the case of py- 
rochlore materials discussed above in the Fig. 6. There- 
fore, as above the long range large interacting moments 
here can be replaced by the interacting tensor fields. It 
was found that the magnetic anisotropy and arrange- 
ments of magnetic moments in BiCoPOs can be con- 
tinuously tuned by temperature and external magnetic 
field [42]. 

However, neither field induced order-disorder nor spin- 
flop transition was found in susceptibility and specific 
heat measurements under external field, although the 
spin-flop transition was seen in the magnetization depen- 
dence on magnetic field measured at T=2 K [12]. The 
Néel temperature, Ty and AF order were found to de- 
crease gradually with the increase in field. But even at 
high fields no extra features other than the Néel transi- 
tion was observed in specific heat [12] until the AF order 
was completely suppressed. Such behavior is consistent 
with the prediction of our holographic model suggest- 
ing that the quite strong and continuous change in mag- 
netic anisotropy with the field may be well described with 
the interacting tensor field embedded in the gravitational 
AdS space, that is with the aid of the AdS/CFT duality. 
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